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ABSTRACT 

Aims. Our aim is to obtain more information about the physical nature of B-type asteroids and 
extend on the previous work by studying their physical properties derived from fitting an asteroid 

in ' 



> 
oo 



thermal model to their NASA's Wide-field Infrared Survey Explorer (WISE) data. We also exam- 
ine the Pallas collisional family, a B-type family with a moderately high albedo in contrast to the 



CO 

o 

CO ' large majority of B-types. 

Methods. We apply a combination of the Near-Earth Asteroid Thermal Model and a model of 

the reflected sunlight to WISE asteroid data in order to derive up to four parameters: effective 

. diameter (£)), the so-called infrared beaming parameter (of), ratio of infrared to visible albedo 

5—i < 

. (ftp ~ PirIPv) an d visible geometric albedo (p v ). 

Results. We obtained the effective diameter, geometric visible albedo, infrared-to-visible albedo 
ratio and beaming parameter for ^ 100 B-types asteroids and plotted the value distributions of 
pv, R P and t] (p v = 0.07 ± 0.03, R,, = 1.0 ± 0.2, and 77 = 1.0 ± 0.1). By combining the IR and 
visible albedos with 2.5 p.m reflectances from the literature we obtained the ratio of reflectances 
at 3.4 and 2.5 pm, from which we found statistically significant indications that the presence of 
a 3-pm absorption band related to water may be commonplace among the B-types. Finally, the 
Pallas collisional family members studied (~ 50 objects) present moderately high values of p v 
(pv = 0.14 ± 0.05), significantly higher than the average albedo of B-types. In addition, this 
family presents the lowest and most homogeneously distributed ^,,-values of our whole sample, 
which shows that this group is clearly different from the rest of B-types, likely because its mem- 



1 



V. Alf-Lagoa et al.: Physical properties of B-type asteroids from WISE data 

bers are pieces likely originating from the same region of (2) Pallas, a particularly high-albedo 
B-type asteroid. 

Key words. Minor planets, asteroids: general - Surveys - Infrared: planetary systems 



1. Introduction 

The study of asteroids is important to gain knowledge about the origin and evolution of our plan- 
etary system. Asteroids are relics of the Solar System's formation and the building blocks of 
the terrestrial planets. Primitive asteroids, i.e. those belonging to the so-called spectroscopic C- 
complex and having in general visible geometric albedo pv ^ 0.1 and featureless, flat visible 
spectra, are particularly relevant in this context. Thought to have formed further away from the 
sun than the other asteroid classes, primitive asteroids have experienced less heating and alter- 
ation processes and have a more pristine composition, potentially preserving crucial information 
about the early Solar System. In addition, primitive asteroids play an important role in current 
exobiological scenarios as they delivered complex organic molecules to the early Earth. This or- 
ganic matter is prerequisite for the synthe sis of pre-biotic biochemical compounds that would sub- 
sequently lead to the emergence of life (lMaurettdl2006l and references therein). For these and 



other reasons, upcoming sampl e return space missions have selected prim itive asteroids as primary 



targets: NASA's OSIRIS-Rex (Campin s et al 



( B aru cci et al.ll2012 



de Leon et al 



2010b; 



2Q1 1[) ; and JAXA's Hayabusa 



Lauret ta et al. 

& 



2010), ESA's Marco Polo-R 



The taxonomic classification of primitive asteroids has been traditionally based on their low 
visible albedo ($0.08-0.1), relatively flat or slightly blue visible spec tra and the weak or no ab- 



Clark etal 



2010). Several pri mitive classes 



1984 



1989). Bus' feature- 



sorption features thereof (for a detailed review, see 
were defined in Tholen's taxonomy, e.g. B, C, F, G, D and P (iTholen 
based classification, independent of the albedo, merged some of these and d efined new primitive 
taxons that were ex tended with minor changes into the Bus-DeMeo taxonomy ( Bus & BinzeJ2002 : 



DeMeo et al.ll2009l) . 

B-types are of particular interest among the primitive asteroids for a number of reasons: (i) 
there is as yet no compelling explanation for their defining feature, i.e. their slightly blue spectral 
slope in the visible range; (ii) B-type asteroids constitute the only primitive class that presents a 
wide range of spectral slopes in the 0. 8-2.5 ^m near-infr ared region, from negative to positive, 
in the 0.8-2.5 fim near-infrared region (Ide Leon et al.ll2012l) ; (iii) the few B-types studied present 
the 3-um absorption feature related to hydrated minerals; (iv) water ice has been detected on the 

j — ' inn i ^ — i 

surface of (24) Themis (ICampins et al.l l2010a; Rivkin & Emery 2010); (v) the majority of aste roids 



that have been observed to display cometary-like activity are B-types (Licandro et al 



2012 



and 



references therein); (vi) the target of NASA OSIRIS-Rex mission, 2006 RQ36, is a B-type asteroid 
B-type asteroids have been widely related to carbonaceous chondrites, composed of carbona- 
ceous minerals and phyllosilicates, in terms of their generally low albedo and broad spectral prop- 



erties (iGaffev et al. 



1989; 



Vilas & Gaffev 



1989; 



Vilas et al 



1994). More recently, 



de Leon et al 



(12012 ) examined visible to near-infrared (VNIR) spectra of a sample of 45 B-types and found that 
the characteristic negative spectral slope in visible wavelengths diverges into a continuum of grad- 
ually varying spectral slopes in the 0.8-2.5 jum (NIR) range, from a monotonic negative (blue) 



http : //www . j spec . j axa . jp/e/ acti vity/hayabusa2 . html 
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de Leon et al 



(2012) classified their spectra into six "average spec- 



slope to a positive (red) slope. 
tra" or "centroids" representative of the whole sample by means of statistical clustering analysis 



(IMarzo et al 



2009) . These cen troids were compared against meteorite spectra from the RELAB 



database dPieters & Hiroill2004l) . The best meteorite analogues found for the six "clusters" were all 
carbonaceous chondrites with a gradual change in their degree of hydration, from aqueously altered 
CM2 chondrites for the reddest cluster, to the heated/thermally metamorphosed CK4 chondrites for 
the bluest one. 



de Leon et al 



(12012b . which are part of 



This work is an extension of the results obtained by 
a ongoing programme devoted to improving our knowledge of B-types. Our aim here is to study 
the physical properties of B-type asteroids that can be derived by fitting a thermal model to their 
NASA's Wide-field Infrared Explorer (WISE) observations, i.e. effectiv e diameter, beaming p aram- 



eter and pm/pv, where pm is the albedo at 3.4-4.6 pm as defined in 



Mainzer et al 



closely follow the methodology of 
described in Sects. [2]and[3] 



Mainzer et al 



(1201 lbl) . We 



(1201 lbl) . though with a number of differences, as 



Mainzer et al 



d201 Id) studied all groups of spectrophotometrically classified asteroids in the 
Tholen, Bus and Bus-DeMeo taxonomies observed by WISE, including B-types. Those authors 
present visible and NIR geometric albedo distributions and median values of B-types and con- 
clude that, in spite of having analogously low albedos, B-, C-, D- and T-type astero ids can be 

(1201 id) point out 



Mainzer et al 



discriminated from their values of NIR reflectance. In particular, 
that B-types have a lower pm/pv ratio than C-types and attribute this to their blue VNIR slopes 
likely extending o ut to 3—4 pm. Our definitio n of B-type asteroid in this work is different: following 



Clark et al 



fcOlOh and 



de Leon et al 



(2012), we consider all objects that have a flat to slightly blue 
spectral slope in the visible range, i.e. any object that has ever been classified as a B-type, includ- 
ing Tholen's F-types and ambiguous designations. This criterion produces a total of 162 asteroids 
classified as B-types. 

We also study the collisional family of (2) Pallas (hereafter PCF and Pallas, respectively). The 
PCF is interesting for the following reasons: (i) it is a B-type family, given t hat Pallas and the 



very few family me mbers that have been taxonomic ally classified are B-t ypes (IGil-Hutton 



2006 



de Leon et alJ uOlO). and the five members studied in 



de Leon et al 



d2012l) were spectrally related 
to carbonaceous chondrites, which establishes their prim itive nature; (i i ) the average values of 



Masiero et al 



d201ll) are roughly 0.15 



geometric albedo of members of the family calculated by 
(see their Fig. 19), significantly greater than expected for primitive bodies (< 0.1), though no 
explicit comment is made by these authors on this intriguing result; (iii) the Near-Earth asteroid 
(3200) Phaethon, an activated asteroid parent of the Geminid meteor shower, likely originated in 
the Pallas family dde Leon et alJuOlOl) : (iv) this family is well isolated in (proper) element space, 
thus the potential identification of interlopers as members is greatly reduced. 

The paper is organised as follows. In Sect. [2] we briefly describe the WISE data set and 
our selection criteria. The thermal modelling of the data is explained in detail in Appendix lAl 
whereas Sect. [3] includes relevant comments on the very few difference s introduced in this work. 



Masiero et al 



(2011). We 



In Appendix IB"! our parameter determinations are compared to those by[ 
present our results in Sect. |4] a discussion of the implications of this work is put forward in Sect. 
[5] and our conclusions are enumerated in Sect. [6] 
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2. Data 



A general introduction to WISE can be found in lWright et al.l (120101) and references therein. Of par- 
ticular interest to Solar System science is the NEOWISE project. This acronym collectively refers 
to two enhancements to the WISE data processing s ystem that were desig ned to allow detection 



Mainzeret al. 2011a). 



and archiving of Solar System objects (for details, see 

WISE used four broad-band filters with approximate i sophotal wavelengt hs at 3.4, 4.6, 12 and 
22 /Km, referred to as Wl, W2, W3 and W4, respectively ( Wright et allboiO ). The WISE All-Sky 
Single Exposure Lib Working Database, published in April 2012 and available via the IRSA/IPAC 
archival includes the corresponding magnitudes and uncertainties in the Vega syste m as well as 



qualit y and contamination and confusion flags that enable us to reject defective data (iCutri et al 
20121) . 

J2OI lblQ) : 



We follow a combination of criteria found in 



Mainzer et al 



Masiero et al. 



and 



Gray et al 



(2011) 



(2012) in order to ensure the reliability of the data. We implement the correction to 
the red and blue calibrator discrepancy in W3 and W4; we use a cone search radius of 0.3" centred 
on the MPC ephemeris of the object in our queries; all artifact flags other than p, P and and 
quality flags other than A, B and C are rejected; we require the modified Julian date to be within 
4 seconds of the time specified by the MPC and split groups of epochs separated more than three 
days (see the end of this section); we ensure that the data is not contaminated by inertial sources 
by removing those points that return a positive match from the WISE Source Catalog within 6"; 
finally, all remaining observations in a given band are rejected if they are fewer than 40% of the 
data in the band with the maximum number of detections. 

On the other hand, we do not use data saturated to any extent. The onset of saturation is reported 
to correspond to magnitudes Mwi < 6, Mw2 < 6, Mw3 < 4, Mw4 < 3 ( Cutri et al]|2012 ). We found 
that enlarging the error bar of partially saturated data to .2 magnitudes (which translates into a 
relative error of 20% in fluxes) renders the corresponding band to play no effective role in the 
thermal model fit by not contributing significantly to the x 1 - 

The application of the above criteria results in a sample of 1 1 1 B-type Main-Belt objects with 
WISE observations usable for our purposes. Some asteroids have been observed by WISE in more 
than one uninterrupted group of epochs with different observation geometries. We also model sep- 
arately such groups of observations if they are more than three days apart (see Appendix [A] and 



Mainzer et al 



201 lbl) . Consequently, we have a larger set of parameter determinations than aster- 



oids in the sample. 



3. Thermal modelling 

The mod e lling of WISE astero id data implemented in this work closely follows 



Mainzer et al 



Masiero et al 



(l2011bEh:l1 

Harrisll 1998b and the IAU phase curve correction to the visible magnitude dBoweil et al.ll 19891) . For 



(201 1) and is based on the Near-Earth Asteroid Thermal model (NEATM, 



the sake of reproducibility, we include a detailed account of our procedure in Appendix [A] and 
enumerate the few differences with respect to Masiero et al. and Mainzer et al. below. 

The number of parameters we can fit for each object depends upon how many and which 
WISE bands are present in its data set. Parameter default values are chosen based on the peak of 



http://irsa.ipac.caltech.edu/Missions/wise.html 
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-0.6 -0.4 



-0.2 

AH 



0.2 0.4 0.6 



Fig. 1. Differences in absolute magnitude values between those used bv lMasiero et al 



(2011), H M , 



and the most updated ones (as of May 2012) used in this paper, H\j. Note that the cases verifying 
AH = do not contribute to this histogram. 



their respective fitted value distributions of Main Belt Asteroids presented in lMasiero et al.l (J201 1 ). 
Whenever there is one single or no thermal band available (W2, W3 or W4) we assume 77 = 1.0; 
R p is is fixed to 1.5 unless we have at least 50% contribution of reflected sunlight in Wl data. The 
last criterion is based on the consistency of our parameter determinations for objects with dou- 
ble detection^. Namely, both groups of observations of asteroids (1076) and (2446) have > 50% 
sunlight and their R p values are consistent within the errorbar. On the other hand, (3579) has non- 
compatible R p determination from Wl data with >70% and ~ 25% reflected sunlight, respectively. 
We thus reject three R p values belonging to asteroids (288), (1493) and (3579). 

These considerations allow asteroid size to be fitted in all cases and, by means of the relation 

lQ -H/5 v2 



py= 1329 [km]—— , (1) 
\ D [km]/ 

the geometric visible albedo can be computed. In contrast, Wl and W2 data are more often rejected 
based on the data requirements (see Sect. |2]i than the purely thermal bands and one will usually be 
able to obtain fewer R p determinations than 77 or indeed D. 

It is important to point out that we do not use physical data previously determined by direct 
measurements such as radar diameters or albedos to constrain our fits, i.e. we limit ourselves to 
using radiometrically derived sizes. This will introduce variations on the parameter determinations 
of some individual asteroids -specially the largest o nes since it is more l ikely that more direct 
measurements have been performed- as compared to 



Masiero et al 



(2011), but should not affect 



the result of statistical analyses if the populations studied consist of a significant number of objects. 
In Appendix [B] we show that, given the same input values of H, our best-fit parameter values are 
consistent within the errorbars, though we find that our R p determinations are systematically lower 
by ~ 10%. On the other hand, the update of ~ 45% of the MPC //-values in our sample does change 
the values of py (see Eq.[TJ and R p . As shown in Fig.Q] the updates tend to be t oward greate r value s 



of H, which will result in lower py and greater values of R p than those of iMasiero et al.1 (120111) . 
For more details, see Appendix IE1 



3 This criterion may not be of general applicability and has only been checked for objects in our B-type 
sample. 
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Fig. 2. Histograms of value distributions obtained for the B-type asteroids observed by WISE: 77 
(top left), p v (bottom left), R„ (top right), and p IR (bottom right). The corresponding histograms 
using the results from 



Masiero et al 



(1201 ll) are overplotted. 



4. Results 

4.1. Value distributions ofrj, p v andR p 

The distributions of 77, py, R p and pm obtained for the B-types are shown in Fig. [2] (the complete 
set of parameters determinations is given in Table [1}. Fo r comparison, we ove rplotted the corre- 
sponding histograms with best-fit parameter values from [Masiero et al.1 (1201 lb (see Appendix [B] 
for a detailed comparison). The parameter median and mean values and standard deviations of this 
work, as well as the number of parameter determinations obtained in each case (AO, are presented 
in Table [2] Note that py and pm are not fitted, but computed. The former is obtained from Eq. Q~| 
with the best-f it value of D as input, whereas pm = R p pv- These results are consistent with pre- 
vious work by 



Mainzer et al. 



(201 lc): if we take the weighted mean of median py- and /? p -values 
corresponding to their Tholen, Bus and Bus-DeMeo B-types and Tholen F-types we obtain the 
same median values. 

We find a 77-value distribution centered at unity, consistent with the average value obtained for 
the whole Main Belt dMasiero et alj|201 11) . The broad and asymmetrical py distribution extends 
to py > 0.1. The R p distribution is also broad, whereas the values of pm are more compactly 
distributed around the mean. 



4.2. Albedo ratio 

Figure [3] shows a plot of R p versus py. Similar plots includ ing all taxonomic clas ses in different 
classification schemes are presented in Figs. 14 and 15 bv Mainzer et al.1 (1201 lch to show how 
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Table 2. Median and mean values and standard deviations of rj, R p , py and put derived for the 
B-type asteroids observed by WISE. 



Parameter 


Median 


Mean 


a 


N 


n 


1.0 


1.0 


0.1 


116 


R P 


1.0 


1.0 


0.2 


88 


Pv 


0.06 


0.07 


0.03 


132 


PlR 


0.06 


0.06 


0.01 


88 



clearly different taxons may be distinguished. Here we concentrate on the py < 0.18 range, with 
all B-types with WISE data for which /^-values could be derived are plotted in black circle s (see 
Table 13; we also include all main belt asteroids taken from Table 1 of iMasiero et al.l (1201 II) (grey 
empty circles). The cloud of points exhibits a characteristic "waning-moon" shape, with no points 
in either the high- py, high-/? p or low-py, low-R p regions of the plot. Mainzer et al. caution that, 
while WISE is essentially unbiased against py, spectroscopic surveys conducted to create the clas- 
sification schemes are inherently biased against small, low-py objects; in addition, the computa- 
tion of R p from WISE data requires sufficient reflected sunlight contribution in bands Wl and W2, 
which will tend to exclude objects with low enough values of py and R p . This could explain the 
lack of points in the lower left part of the plot. However, these biases cannot be solely responsible 
for this characteristic shape since other taxonomic classes with higher values of pv and R p also 
cluster similarly. Furthermore, if we plot all main belt objects irrespective of whether they have 
a taxonomic classification or not, distinct clouds of points with the same shape become apparent. 
Thus, we emphasize that because of their characteristic py- and /?/R-value distributions, there are 
no high-/?,, objects among the high-py B-type asteroids. 



2.5 
2 

1.5 
a* 1 

0.5 


-0.5 









All main belt asteroids 

B-tvpes • 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 

Pv 



Fig. 3. Albedo ratio versus visible geometric albedo. The 20% errorbars in py and R p are not 
shown to facilitate visualisation. B-types observed by W ISE are plotted in black circles; all main 
belt objects featured in Table 1 of 



Masiero et al. 



(1201 II) are in empty grey circles. The horizontal 
line is an artificial feature corresponding to objects with the default fixed value of R p -\.5. 



We have also analysed the /^-values of the sample of 45 B-types studied by 



de Leon et al 



d2012l) separately. The spectra of these asteroids were classified into six "average spectra" or "cen- 
troids" referred to as Gl, G2, ... G6 (see Sect. [Tj. These show a progressive decrease in spectral 
gradient in the NIR interval (0.8-2.5 pm), ranging from a positive (red) slope for Gl to a negative 
(blue) slope for G6. In Fig. [4] we plot R p versus p y labelling the objects in the different centroids 
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G4 (0.08 ± 0.03) 
G5 (0.08 ± 0.02) 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 

Pv 



de Leon et al 



(2012) distinguishing the clusters to which they 



Fig. 4. R p vs. pv for the asteroids in 
belong. Within parentheses, the average py of each cluster is shown. On average, pv increases from 
Gl to G3, keeping the same value from G3 to G5. 



1.4 
1.3 
1.2 
aT 1-1 

tot) j 

2 1 
< 0.9 
0.8 
0.7 
0.6 



G1 



G2 G3 G4 

Cluster 



G5 



de Leon et al 



(2012). In terms of spectral 



Fig. 5. Average values of R p for each cluster defined by < 
slope up to 2.5 pm, cluster G5 is the bluest, whereas Gl is the reddest. Note: the errorbar is the 
standard error of the mean. 



Gl, . . . G5. Note that cluster G6 is not included since it is composed only of one member, (3200) 
Phaethon, which did not have enough WISE observations to perform a reliable fit. The average 
values of py for each cluster increases from Gl to G3, retaining the same value from G3 to G5. 
This might suggest an inverse correlation between the cluster NIR slope and py, though the small 
number of objects per cluster with WISE observations (2 objects in Gl, 12 in G2, 6 in G3, 5 in G4, 
and 3 in G5) prevents us from stablishing a firm conclusion. 

In Fig. [21 the average /? p -value is plotted for the different clusters. This figure suggests a corre- 
lation between the average R p and the NIR slope of the relative reflectance of the clusters: objects 
with higher R p belong (o n average) to cluste rs with higher NIR spectral slope in the ~ 1.0-2.5 pm 
range (see also Fig. 5 of Ide Leon et al I2012L This indicates that the refle ctivity at 3.4 pm tends to 



continue the trend observed at shorter IR wavelength s, as hypothesized by 



Mainzer et al 



The Wl bandpass spans the 2.8 to 3.8 pm range (Wrig ht et al. 



(1201 id) . 



2010); therefore, the values of R,- 



may also be diagnostic of the presence of the 3-pm absorption feature attributed to hydrated miner- 



als or water-ice detected on many asteroids ( Ri vkin 



2010b; 



Rivkin & Emery 



2010; 



Licandro et al. 



i et 



2011). 



al. 



2000; 



Gaffev et al 



2002; 



Mainzer et alj d201 Id) ruled out the possi 



Campins et al 



bility of detecting the hydration band from WISE data base d on the fact that t he average R p for 



a sample of 7 M-types with positive detections of the band dRivkin et al 



2000) cannot be distin- 



guished from that corresponding to other 33 M-types. However, this test might not be meaningful 
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Fig. 6. Assembled spectra of the control objects chosen to study the reflectivity at 3.4 yum derived 
from WISE data as diagnostic of the presence of the 3-pm feature. The dashed lines approximately 
enclose the wavelength integration range of the Wl filter. 



given that the /? ; ,-values have errorbars at least larger than the characteristic depth of the absorption 
feature. Below we provide evidence that the Wl may be sensitive to the 3-yum feature. 

Our first step was to assemble VNIR spectra up to ~ 3.6 pm of a list of nine C-complex control 
asteroids. All of these asteroids are primitive, including some B-types, such as (2) Pallas or (45) 
Eugenia, and some of them show a distinct absorption feature. The assembled spectra are plotte d 
along with t he best-fit value of R D in Fig. [6] We collecte d or digitised dat a from 



SMASS-II (Bus & Binzel 2002), the 52-Color Survey ( Bell et al 



de Leon et al 



2005), 



Hiroi et 



Rivkin et al 



al 



1996), 



(12003b and 



(2012). Taking into account the errors associated with superimposing spectra ob- 
tained at different epochs of observation with different equipment and the uncertainties in the spec- 
tra and in the value of the albedo ratio, the latter is not expected to exactly match the value of the 
spectra at 3.4 pm. Owing to these deviations, seen in Fig. [6] one cannot confirm nor rule out the 
presence of the absorption feature based on the value of R p alone. 

Next, we combi ned the values of R„ with the observed relative reflectances at 2.5 pm available 

(2003) use the parameter 1 - Rx/Ri.s as a rough measure of band 



from other datasets. 



Rivkin et al 



depth. Because the band minima are usually near 3.0 pm and given that Wl results from an average 
over 2.8 to 3.8 pm, the parameter b = 1 - R p /R2.5 is not to be taken as a measure of band depth but 
as a helpful parameter to quantitatively compare the values of relative reflectances at 2.5 and 3.4 
pm. If a deep absorption band is present, R p is in general expected to be S R2.5 and hence b ^ 0, as 
is the case of (19) Fortuna. There is one general case for which this interpretation would be wrong: 
if the NIR slope is negative up to 3.4 pm, we would have b ^ even if no band was present, though 
blue NIR slopes such as Pallas' have only been measured in a small percentage of cases (e.g. 2 out 
of 45 in the sample of 



de Leon et al. 2012). 
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Fig. 7. Plot of b — 1 - R p /R2.5 versus pv for a set of primitive control objects. The values of b 
obtained from the combination of our R p and the R2.5 taken from the assembly of spectra published 
by other authors (cf. figure |6) and the shape of the 3 pm absorption band are consistent except for 
(13) Egeria. 



A plot of b vs. pv is shown in Fig. [7] with an estimated b errorbar of 0.2 to account for the 
large uncertainties in R p and the assembly of the spectra. We see that those asteroids with a weak 
or non-existent absorption band tend to have b S= 0, whereas b ^ for those with a higher-contrast 
feature. Only in the case of (13) Egeria the value of b is clearly inconsistent with the absortion 
band observed. Indeed, the large errorbar in the value of b still prevents irrefutable detection of the 
band on a case-by-case basis, but from Fig. [7] the correlation between the sign of b and the spectra 
appears robust, so that systematically obtaining b ^ for a given population may be statistically 
significant. 

In Fig. [8] (left panel) we show a fe-value histogram of those B type asteroids for which the 
value of b could be determined. R2.5 was taken from the B-type spectra presented in Ide Leon et al 



d2.012l) . after normalising to unity at 0.55 pm. In order to test for statistical significance, we resort 
to the Kolmogorov-Smirnov (KS) test, which enables one to reject the null hypothesis that a given 

set of unbinned values is compatible with having been drawn from a given distribution function 

I II 1 

(see e.g. IPress et al. 1986). As the null hypothesis we take a gaussian distribution with zero mean. 
This choice is based on the mean value of b that we would expect considering that interpolating 
between R0.55 = 1 and R p = 1 gives R2.5 = 1-0 => b = (see Fig. |9]». The KS test amply rules out 
the null hypothesis that the B-types ^-values are drawn from a gaussian distribution of zero mean 
regardless of its width (<x). 

To demonstrate the robustnest of this result even further, we carried out the same procedure for 
a list of asteroids belonging to the S-complex, including all Bus-DeMeo pure S-types for which 
we found the value of R2.5 was available from the literature and for which R p could be computed, 
plus enough randomly selected S subtypes to get the same number of /b-values we derived for the 
B-types (see Table©. S-type asteroids are "anhydrous" and have positive spectral slopes up to 2.5 
pm so that in the absence of a 3-pm absorption feature one would expect a negative value of b if 
the spectral slope maintains its trend up to 3.4 pm. Taking into account that R p =s 1.67 for this 
sample of S-types, from the interpolated value R2.5 — 1-46 (see Table [3}, the expected mean value 
of b would be b — -0.14 (see Fig.[9]l. The distribution of ^-values obtained for the S-types, shown 
in the right panel of Fig. [8] presents a negative mean value. In this case, the KS test does not rule 
out the null hypothesis that the /b-values are drawn from a Gaussian distribution centered at -0.14 
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Fig. 8. Normalised histograms showing the distribution of ^-values. Left panel: B-types; the fact 
that b > in the majority of cases suggests that a large percentage of these asteroids have absorption 
features in the ~ 3-//m region. Right panel: S-types; the opposite conclusion is reached, consistent 
with the anhydrous nature of S-types. 
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Fig. 9. Schematic diagram illustrating the interpolation of R2.5 (empty points) from the average 
values of R p (filled points) for the B and S asteroids studied. 

with cr = 0.2 with p-value > 0.9. This rules out the possibility that a systematic error is causing the 
^-values of the B-type asteroids to be > 0. 



Table 4. Median and standard deviation of fc-values for B-types and S-types 



Group 


Median 


<r 


N 


B-types 


0.2 


0.2 


34 


S-types 


-0.2 


0.2 


34 



In the prec eeding subsection w e showed that our values of R p are systematically ~10% lower 



than those of 



Masiero et al 



d201 II) . We carried out the same procedure enlarging our R p values 
10% and the conclusions still hold. 

To conclude, we find that the majority of B-types with computed ^-values verify b > and 
that very few present a clearly negative value of b, which means that even for thos e B-type spec 
tra w ith a positive slope in the 2.5 fim region (approximately half of the objects in 



de Leon et al 



2012) there is a reduction in the reflectivity around 3.4 fim. As discussed above, for asteroids of the 



C-complex, the 3-fim absorption feature has been attributed to hydrated minerals or wat er-ice. The 



presence of goethite has been also proposed as an alternative explanation for this band (Beck et al 
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2011). Nonetheless, while other closely related minerals have been found in both meteorite and 
asteroid spectra, extrarrestrial goethite has never been identified within the meteorite inventory, 
so the possibility that putative goethite-containing asteroids never found a dynamical collisional 
pathway t o Earth is less likely than the sim pler interpretation: goethite is not present in asteroidal 



surfaces (IJewitt & Guilbert-Lepoutre II201 2I. Therefore, from the distribution of ^-values, we con- 
clude that most asteroids in this sample (which constitutes 5 40% of the B-type population with 
computed R p ) present this absorption and that "water" (be it bound or free) may be common among 
the B-type asteroids. 

4.3. The Pallas collisional family 



2012b- WISE has observed 46 of the 



We use the most up-to-date Pallas family list by (Nesvornv 
objects in this list. Histograms of beaming parameter and albedo determinations are shown in Fig. 
[10](the complete set of parameters is shown in Table|4). Given that the R p value could only be fitted 
for seven objects, we do not include histograms for R p and In Table|6]we present the mean 
values of rj, p v and R p along with their corresponding standard deviations and contributing number 
of determinations. 
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Fig. 10. Beaming parameter (left) and geometric albedo (right) distributions of the values derived 
for the PCF members observed by WISE. The albedo distributions are normalised, i.e. divided by 
their respective total number of counts. 



Table 6. Mean values and standard deviations of rj, R p and p v derived for the members of the Pallas collisional 
family observed by WISE. 



Parameter 


Mean Value 


<r 


N 


1 


1.0 


0.2 


46 


Rp 


0.5 


0.1 


7 


Pv 


0.14 


0.05 


50 



The averag e rj value of the PC F is consistent with that of the B-type population and that of 



the Main Belt (Masiero et al. 



201 1). However, the PCF has a moderately high albedo of ~ 14%, 



significantly higher than the rest of the B-types and indeed higher than the value expected for prim- 
itive asteroids, whereas their average R p is lower than the average value for the B-types (cf. Table 
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Fig. 11. Left panel: PCF members are on average brighter in the visible and bluer in the 3.4 /mi 
region than the B-types sample (cf. Fig. O. The 20% errorbars are not plotted for clarit y. Right 
panel: the average R p of the PCF continues the decreasing trend of the 
(cf.Fig.H). 



de Leon et al 



(2012) clusters 



13. The seven ./^-values computed are quite homogeneous compared to the distribution observed 



for the rest of the B-type population. As we saw in Sect 



consistent within the errorbars with those of (Masiero et al 



4.11 our parameter determinations are 



201 II) . so we do not perform the same 



comparison for the PCF members. Masiero et al. show a histogram of log pv with peak between 
log pv = (-1) - (-0.8), i.e. pv =0.10-0.16, with which our results are consistent (they provide 
no mean value to compare with). On the other hand, though they provide the best-fit parameters in 
their Table 1, Masiero et al. do not discuss the R„ values of the PCF. 



The average pv and R p values continue the trend observed for the clusters of 



de Leon et al 



(12012 ) in Fig. [4] and Fig. [5] they have higher visible albedos on average than the G5 cluster and are 
bluer at 3.4 pm than the bluest cluster, G4 (note that the average R p value of cluster G5 is obtained 
from only three object s, which could exp lain its deviation from the trend). Members of the Pallas- 



de Leon et al 



like gr oup of objects of lClark et al.l (l2010) were concentrated in clusters G4 and G5 in 
(120121) (see Sect. [TJ. In this work, the list of Pallas family members observed by WISE includes 
yet more aste r oids. The four additional objects with R p determinations that were not present in 



de Leon et al 



(120121) are also located in the lower right region of Fig. [4] This property would be 



consistent with the characteristics of an extrapolated G 6 centroid (see the right panel of Fig. ITTb 
Interestingly, the G6 centroid of 



de Leon et al 



(2012) contained asteroid (3200) Phaethon alone. 
Th erefore, the connectio n of the NEA (3200) Phaethon to the Pallas Collisional Family stablished 
by Ide Leon et al.l (120101) based on spectroscopical and dynamical arguments is also supported by 



the values of relative reflectances at 3.4 pm. 



5. Discussion 



In Sect. 14.21 we present statistically significant indications that water may be common within our 
sample, and in Sect. I4.3l we point out how the PCF have higher geometric visible albedos and more 
homogeneously distributed /? p -values than the rest of B-types. The purpose of this section is to put 
these results in the context of other studies in the literature. 

In the context of the geop hysic al models of the Themis and Pallas parent bodies by 
Castillo-Rogez & Schmidll (120 101) and ISchmidt & Castillo-Rogez (2012) , the detections of water 



ice and organics on the surface of (24) Themis (IRivkin & EmervlfeOlOt ICampins et al 



2010bJ) and 
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(65) Cybele (iLicandro et al.l l201 1) have been invoked as provocative indications that water has 
played an important role in the accretion and evolution of these asteroids' parent bodies in the mid- 
outer asteroid belt. Finding indications that a 3-jt/m is frequently present in B-types spectra adds 
further support to this conceptual framework, in which the gradual differences in the observed NIR 
spectral slopes of B-types might also be explained as resulting from different subsequent evolution 
an d processing of their parent bodies . 



Schmidt & Castillo-Rogea (120121) discuss that (1) Ceres, being the largest body of the asteroid 



belt and plausibly a water-rich asteroid, mostly preserved its integrity, whereas the Themis parent 
body was catastrophically disrupted. Pallas would be an intermediate case, still intact but showing 
evidence for heating, internal evolution, impa ct and loss of water. I n accordance with this scenario, 
the Themis family members included in the 



de Leon et al 



(2012) sample were distributed more 
or less homogeneously among the representative centroids of the complete sample, spanning from 
neutral to red slopes and matching several meteorite analogues, whereas the PCF members were 
concentrated in the bluest centroids. As our results show, these trends are also verified longward 
of 2.5 /mi: the Themis family members show more heterogeneous R p values whereas the Pallas 
family members have distinctly higher albedos and lower, more homogeneous values of R p than 
the rest of B-types. This homogeneity is consistent with the collisional family being the result of 
a non-disruptive cratering event on Pallas, which has a significantly higher albedo than the rest of 

large B-types. 

The geophysical models by ISchmidt & Castillo-Rogea (1201 21) favour a water-rich past for 
Pallas. Its near spherical shape may be explained as a consequence of early melting of a substantial 
initial icy component that is subsequently removed by a combination of thermal and impacts pro- 
cesses. The water-loss processes that took place on the surface of Pallas may have taken place on 
other active bodies such as main belt c omets or indeed (3200) Phaethon, linked to the Pallas col- 
lisional family bv lde Leon et alJ (120101) , and this activity has already been proposed to explain the 



characteristic NIR blue spectral slope (ISchmidt & Castillo-Rogez 



2012, and references therein). 



According to this model, water-rich materials are also expected to be important components within 
some members of Pallas' family. The fact th at the PCF members in this study are found to match 
the properties of the least hydrated clusters of c 



de Leon et al 



(2012) would be inconsistent with this 



view, but caution must be exercised when associating NIR slopes of meteorite analogues to hydra- 
tion. For instance, while having suffered intense thermal metamorphism, CK4 chondrites (the best 
matches for Pallas) still show an absorption feature at 2.9 /mi, albeit less prominent than e.g. CM 
chondrites. 



From the dynamical standpoint, the simulations of lWalsh et alJ d201 11) also provide a congruous 
context for our results. Their model suggests that C-complex asteroids were formed in the giant- 
planet-forming region and that, for every C-type planetesimal from beyond 8 AU that would later 
be located in the outer Main Belt, > 10 would have ended up in the region where terrestrial planets 
formed. Assuming that the composition of these objects is 10% water by mass, this may account 
for the minimal mass required to bring the current amount of water to our planet by a factor of 
6-22. 

On the other hand, the question of the moderately high values of pv obtained for the PCF 
remains unanswered. Finding members of the Pallas and Themis collisional families distributed 



differently among the 



de Leon et al 



(2012) clusters and in Fig.QT|(as discussed above) also leads 
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us to ponder the possibility that a sequential or progressive physical process could explain the 
differences in the value distributions of py and R p . Unfortunately, there are many possible mecha- 
nisms underlying the observed NIR spectral variability of primitive asteroids (e.g. regolith particle 
sizes, space weathering as a function of asteroid famili y age and composition, thermal processing) 
and these are difficult to disentangle (Zi ffer et al 



6. Conclusions 



We have derived values of D, 77, and R p of 1 1 1 B-type asteroids by means of therma l modelling of 
WISE data and updated H values (most of which have increased since the work of 



Masiero et al 



2011 



as shown in Fig.[TJ. Additionally, from H and D we have computed the corresponding values 
of py (Tabled. Our results are in a greement within the errorbars of the model with those previously 
published by 



Masiero et al. 



(1201 11) . However, we obtain a systematic trend of approximately -10% 
discrepancies in the determinations of R p that we cannot explain but is most probably attributable to 
small differences in the tabulated solar flux data necessary to estimate the reflected light component 
at 3.4 pm. This work led to the following conclusions: 

1. We derived the distribution of 77, py and Rp fot the B-type asteroids (see Fig. |2]i and obtained 
the following mean values: fj = 1.0 + 0.1, py = 0.07 ± 0.03, and R p = 1.0 + 0.3. 

2. There are no high-py, high-7? p B-type asteroids (py > 0.10, R p > 1.0; see Fig. [3]). The average 



R p -values of the centroids of 



de Leon et al 



(12.012b clearly decrease from Gl to G5, which 



implies a relationship between the IR slope of the asteroid spectra up to 2.5 pm (Figs.[4]and[5]). 

3. We computed ^-values for a set of B-type asteroids and S-complex asteroids, which present 
b > and b < 0, respectively (Fig. [8j. While the latter result is consistent with objects of the 
S-complex being anhydrous, the former indicates that the majority of B-type asteroids ^-values 
in this study are consistent with the presence of a 3-^m absorption feature usually attributed to 
hydrated minerals or water-ice; therefore, water must have played a key role in the evolution of 
a large fraction of the B-types, supporting recent works discussed in Sect. [5] 

4. We have also studied the Pallas collisional family. On the one hand, the average albedo {py = 
0.14 + 0.05) of this familiy is significantly higher than the average albedo of B-types {py = 
0.07 ± 0.03) and moderately high compared to what is traditionally considered to be the albedo 
of primitive asteroids (< 0.1). On the other hand, the albedo ratio values of the PCF members 
are very low and homogeneous {R p = 0.5 ± 0.1). These results clearly show the fundamental 
differences between the Pallas Collisional Family and the rest of B-types. 

In addition , the connection of the NEA (3200) Phaethon to the Pallas Collisional Family sta- 
blished by Ide Leon et alJ (12010 ) based on spectroscopical and dynamical arguments is also 



supported by the values of relative reflectances at 3.4 pm. 



5. Our re sult s support the scenarios by the ge ophysical models by Castillo-Roaez & Sc hmidt 



(2010) and Sch midt & Casfillo-Rogezl(l2012|) and the simulations of lWalsh et al 



2011b . which 



suggest that water played an important role in the origin of primitive asteroid parent bodies of 
the mid-outer belt. 
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V. Alf-Lagoa et al.: Physical properties of B-type asteroids from WISE data 

Table 1. Best-fitting values of physical parameters determined for the B-types with WISE observations. 
Negative values of tj and or R p indicate that the parameter was not free but fixed to the corresponding positive 
value. Errorbars shown are minimum estimates and correspond to 10% relative error for D and 20% for tj, p v 
and R p . 
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Pv 






R P 


Wl 


W2 


W3 


W4 


R2.5 


00002 


4.13 


0. 


11 


669 ± 


67 


0.09 


+ 


0.02 


-1.0 


± 0.2 


0.8 


±0.2 


3 


3 








0.8 


00024 


7.08 


0. 


19 


187 ± 


19 


0.07 


+ 


0.01 


-1.0 


± 0.2 


1.0 


±0.2 


5 


5 








1.1 


00045 


7.46 


0.07 


240 ± 


24 


0.03 


+ 


0.01 


-1.0 


± 0.2 


1.6 


±0.3 


9 


10 








1.3 


00047 


7.84 


0. 


16 


118 ± 


12 


0.09 


+ 


0.02 


-1.0 


± 0.2 


0.8 


±0.2 


9 


9 










00052 


6.31 


0. 


18 


396 ± 


40 


0.03 


+ 


0.01 


-1.0 


± 0.2 


1.4 


±0.3 


4 


4 








1.3 


00059 


7.93 


0. 


15 


169 ± 


17 


0.04 


+ 


0.01 


-1.0 


± 0.2 


1.4 


±0.3 


11 


11 








1.3 


00085 


7.61 


0. 


15 


148 ± 


15 


0.07 


+ 


0.01 


-1.0 


+ 0.2 


1.1 


±0.2 


6 


6 








1.1 


00141 


8.40 


0. 


15 


129 ± 


13 


0.05 


+ 


0.01 


-1.0 


± 0.2 


1.0 


±0.2 


9 


9 










00142 


10.27 


0. 


15 


58 ± 


6 


0.04 


+ 


0.01 


1.0 


± 0.2 


1.1 


±0.2 


8 


8 





8 


1.1 


00225 


8.72 


0. 


15 


107 ± 


11 


0.05 


+ 


0.01 


1.1 


± 0.2 


1.1 


± 0.2 


14 


15 


15 


14 




00229 


9.13 


0. 


15 


110 ± 


11 


0.03 


+ 


0.01 


0.8 


± 0.2 


1.5 


±0.3 


10 


10 





10 




00241 


7.58 


0. 


15 


189 ± 


19 


0.05 


+ 


0.01 


-1.0 


± 0.2 


1.3 


±0.3 


7 


7 










00241 


7.58 


0. 


15 


198 ± 


20 


0.04 


+ 


0.01 


-1.0 


± 0.2 


1.3 


±0.3 


7 


7 










00268 


8.28 


0. 


15 


142 ± 


14 


0.04 


+ 


0.01 


-1.0 


± 0.2 


1.4 


±0.3 


9 


9 










00282 


10.91 


0. 


15 


41 ± 


4 


0.05 


+ 


0.01 


1.0 


± 0.2 


-1.5 


± 0.3 


10 


10 





10 




00314 


9.80 


0. 


15 


64 ± 


6 


0.05 


+ 


0.01 


1.0 


± 0.2 


1.0 


±0.2 


14 


14 





14 




00335 


8.96 


0. 


15 


88 ± 


9 


0.06 


+ 


0.01 


1.1 


± 0.2 


0.9 


±0.2 


4 


4 





4 


1.1 


00357 


8.72 


0. 


15 


105 ± 


10 


0.05 


+ 


0.01 


1.0 


± 0.2 


0.8 


±0.2 


8 


8 





8 


1.1 


00372 


7.50 


0. 


15 


180 ± 


18 


0.05 


+ 


0.01 


-1.0 


± 0.2 


1.3 


±0.3 


12 


12 










00379 


8.87 


0.15 


88 ± 


9 


0.06 


+ 


0.01 


1.0 


± 0.2 


1.0 


±0.2 


6 


6 





6 


1.2 


00383 


9.91 


0. 


15 


44 ± 


4 


0.10 


+ 


0.02 


1.3 


± 0.2 


0.7 


± 0.1 


12 


13 


13 


13 


0.9 


00400 


10.50 


0. 


15 


39 ± 


4 


0.07 


+ 


0.01 


1.2 


± 0.2 


0.8 


±0.2 


1 1 


1 1 


11 


10 




00404 


9.01 


0. 


15 


98 ± 


10 


0.05 


+ 


0.01 


1.0 


± 0.2 


0.9 


±0.2 


8 


8 





8 




00426 


8.42 


0. 


15 


117 ± 


12 


0.06 


+ 


0.01 


0.9 


± 0.2 


1.1 


±0.2 


9 


9 





6 


0.8 


00431 


8.72 


0. 


15 


103 ± 


10 


0.05 


+ 


0.01 


0.9 


± 0.2 


1.1 


±0.2 


13 


12 





12 


1.2 


00461 


10.48 


0. 


15 


46 ± 


5 


0.05 


+ 


0.01 


1.0 


± 0.2 


1.1 


±0.2 


26 


27 


27 


27 




00464 


9.52 


0. 


15 


84 ± 


8 


0.04 


+ 


0.01 


1.0 


± 0.2 


1.1 


±0.2 


4 


4 





4 




00464 


9.52 


0. 


15 


82 ± 


8 


0.04 


+ 


0.01 


0.9 


± 0.2 


1.2 


± 0.2 


26 


26 





26 




00468 


9.83 


0. 


15 


66 + 


7 


0.05 


+ 


0.01 


1.0 


± 0.2 


1.2 


±0.2 


13 


13 





13 




00526 


10.17 


0. 


15 


48 ± 


5 


0.06 


+ 


0.01 


1.1 


± 0.2 


0.9 


±0.2 


13 


12 


17 


17 




00531 


12.00 


0. 


15 


16 ± 


2 


0.10 


+ 


0.02 


0.9 


± 0.2 


-1.5 


±0.3 





8 


12 


12 




00541 


10.10 


0. 


15 


57 ± 


6 


0.05 


+ 


0.01 


1.0 


± 0.2 


1.0 


±0.2 


13 


11 





13 




00555 


10.70 


0. 


15 


33 ± 


3 


0.09 


+ 


0.02 


1.0 


± 0.2 


0.7 


± 0.1 


12 


12 


12 


12 




00560 


10.90 


0. 


15 


36 ± 


4 


0.06 


+ 


0.01 


1.0 


± 0.2 


0.8 


±0.2 


14 


14 


9 


14 




00567 


9.16 


0.15 


91 ± 


9 


0.05 


+ 


0.01 


1.0 


±0.2 


1.4 


±0.3 


8 


8 





7 




00567 


9.16 


0. 


15 


82 ± 


8 


0.06 


+ 


0.01 


1.0 


±0.2 


1.1 


±0.2 


8 


8 





8 




00635 


9.01 


0. 


15 


97 ± 


10 


0.05 


+ 


0.01 


0.9 


± 0.2 


0.9 


±0.2 


11 


11 





8 




00702 


7.25 


0. 


15 


202 ± 


20 


0.05 


+ 


0.01 


-1.0 


± 0.2 


1.2 


± 0.2 


10 


10 










00704 


5.94 


-0.02 


361 ± 


36 


0.06 


+ 


0.01 


-1.0 


± 0.2 


1.1 


± 0.2 


6 


7 










00704 


5.94 


-0.02 


351 ± 


35 


0.06 


+ 


0.01 


-1.0 


± 0.2 


1.1 


± 0.2 


8 


8 











18 



V. Alf-Lagoa et al.: Physical properties of B-type asteroids from WISE data 
Table 1. continued. 
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V. Alf-Lagoa et al.: Physical properties of B-type asteroids from WISE data 
Table 1. continued. 
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11 ± 1 


0.06 ± 0.01 


1.0 ± 0.2 


-1.5 ± 0.3 








15 


15 




03036 


10.30 


0. 


15 


47 ± 5 


0.06 ± 0.01 


1.0 ± 0.2 


1.0 ± 0.2 


9 


10 


10 


10 


1.3 


03036 


10.30 


0. 


15 


51 ± 5 


0.05 ± 0.01 


1.0 ± 0.2 


1.0 ± 0.2 


8 


8 


8 


8 


1.3 


03074 


13.60 


0. 


15 


10 ± 1 


0.06 ± 0.01 


1.2 ± 0.2 


-1.5 ± 0.3 





5 


9 


9 




03139 


10.70 


0. 


15 


41 ± 4 


0.06 ± 0.01 


1.1 ± 0.2 


0.9 ± 0.2 


10 


12 


11 


10 




03162 


11.30 


0. 


15 


31 ± 3 


0.05 ± 0.01 


1.1 ±0.2 


-1.5 ± 0.2 


13 


8 


16 


16 


1.3 


03162 


11.30 


0. 


15 


33 ± 3 


0.05 ± 0.01 


1.1 ±0.2 


0.9 ± 0.2 


10 


10 


10 


11 


1.3 


03204 


12.20 


0. 


15 


21 ± 2 


0.05 ± 0.01 


1.2 ± 0.2 


-1.5 ± 0.3 








7 


7 




03204 


12.20 


0. 


15 


20 ± 2 


0.06 ± 0.01 


1.3 ± 0.2 


-1.5 ± 0.3 








7 


7 




03566 


12.90 


0. 


15 


14 ± 1 


0.07 ± 0.01 


0.9 ± 0.2 


-1.5 ± 0.3 





5 


8 


8 




03579 


13.60 


0. 


15 


7 ± 1 


0.12 ± 0.02 


0.9 ± 0.2 


0.4 ±0.1 


10 


14 


13 


13 


1.0 


03579 


13.60 


0. 


15 


7 ± 1 


0.13 ± 0.03 


1.0 ± 0.2 


-1.5 ± 0.3 








9 


9 


1.0 


03581 


12.10 


0. 


15 


14 ± 1 


0.13 ± 0.03 


0.9 ± 0.2 


-1.5 ± 0.3 








13 


13 




03627 


13.50 


0. 


15 


11 ± 1 


0.06 ± 0.01 


1.0 ± 0.2 


-1.5 ± 0.3 





10 


10 


10 




03647 


11.50 


0. 


15 


28 + 3 


0.06 ± 0.01 


1.0 ± 0.2 


0.8 ± 0.2 


11 


12 


12 


12 




04100 


11.50 


0. 


15 


17 ± 2 


0.15 ± 0.03 


1.0 ± 0.2 


0.9 ± 0.2 


5 





1 1 


11 


1.1 


04396 


13.60 


0. 


15 


5 ± 1 


0.24 ± 0.05 


1.1 ±0.2 


-1.5 ± 0.3 








8 


8 




04484 


12.30 


0. 


15 


16 ± 2 


0.08 ± 0.02 


1.1 ±0.2 


-1.5 ± 0.3 





6 


9 


9 




04837 


11.60 


0. 


15 


28 ± 3 


0.05 ± 0.01 


0.9 ± 0.2 


-1.5 ± 0.3 








10 


10 




04955 


11.70 


0. 


15 


22 ± 2 


0.08 ± 0.02 


1.2 ± 0.2 


-1.5 ± 0.3 








12 


12 




04997 


12.70 


0. 


15 


10 ± 1 


0.16 ± 0.03 


1.1 ±0.2 


-1.5 ± 0.3 








10 


10 




05057 


12.30 


0. 


15 


18 ± 2 


0.06 ± 0.01 


1.0 ± 0.2 


-1.5 ± 0.3 








8 


8 




05133 


11.90 


0. 


15 


24 ± 2 


0.05 ± 0.01 


1.1 ± 0.2 


-1.5 ± 0.3 





8 


12 


12 




05222 


11.30 


0. 


15 


iy ± z 


0.15 ± 0.03 


l.t) ± 0.2 


U.4 ±0.1 


12 


13 


13 


13 




05234 


11.90 


0. 


15 


15 ± 1 


0.15 + 0.03 


1.0 ± 0.2 


-1.5 ±0.3 








6 


6 




05234 


11.90 


0. 


15 


14 ±2 


0.13 ± 0.03 


0.9 ± 0.2 


-1.5 ± 0.3 








10 


10 




05330 


12.10 


0. 


15 


15 ±2 


0.11 ± 0.02 


0.9 ± 0.2 


0.6 ±0.1 


9 


16 


16 


15 




05870 


13.20 


0. 


15 


9± 1 


0.11+0.02 


1.0 ±0.2 


-1.5 ±0.3 








5 


4 




06297 


12.30 


0. 


15 


18 ±2 


0.07 ± 0.01 


1.1 ±0.2 


-1.5 ±0.3 








8 


8 
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Table 1. continued. 



Designation 


H 


G 


D [km] 


Pv 


i] 


R P 


Wl 


W2 


W3 


W4 R 25 


08518 


12.80 


0.15 


14+ 1 


0.07 ± 0.01 


1.1+0.2 


-1.5 + 0.3 








5 


5 


08519 


13.60 


0.15 


7 + 1 


0.12 + 0.02 


1.2 + 0.2 


-1.5 + 0.3 








11 


11 


08906 


12.80 


0.15 


14+ 1 


0.07 ± 0.01 


0.9 ± 0.2 


-1.5 + 0.3 








6 


6 


09219 


11.90 


0.15 


20 ±2 


0.07 ± 0.01 


1.3 + 0.2 


-1.5 + 0.3 








8 


8 
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Table 3. Best-fitting values of physical parameters determined for the S-types with sufficient Wl WISE obser- 
vations and published 2.5 pm reflectances. Negative values of rj and or R p indicate that the parameter was not 
free but fixed to the corresponding positive value. Errorbars shown are minimum estimates and correspond to 
10% relative error for D and 20% for r], p v and R p . 



Designation 


H 


G 


D [km] 


Pv 




R P 


Wl 


W2 


W3 


W4 


^25 


00009 


6.280 


0.170 


231 ± 23 


0.10 ± 0.02 


-1.0 ± 0.2 


2.1 ± 0.4 


8 


8 








1.5 


00009 


6.280 


0.170 


188 ± 19 


0.15 ± 0.03 


-1.0 ± 0.2 


1.9 ± 0.4 


13 


13 








1.5 


00011 


6.550 


0.150 


171 ± 17 


0.14 ± 0.03 


-1.0 ± 0.2 


1.6 ± 0.3 


8 


8 








1.4 


00017 


7.760 


0.150 


93 ± 9 


0.16 ± 0.03 


0.9 ± 0.2 


1.8 ± 0.4 


10 


10 





10 


1.5 


00026 


7.400 


0.150 


86 ± 9 


0.26 ± 0.05 


1.1 ± 0.2 


1.3 ± 0.3 


11 


12 





12 


1.4 


00029 


5.850 


0.200 


194 ± 19 


0.21 ± 0.04 


1.0 ± 0.2 


1.4 ± 0.3 


10 


10 








1.3 


00030 


7.570 


0.150 


88 ± 9 


0.21 ± 0.04 


-1.0 ± 0.2 


1.8 ± 0.4 


11 


11 








1.5 


00032 


7.560 


0.150 


78 ± 8 


0.27 + 0.05 


-1.0 ± 0.2 


1.4 ± 0.3 


7 


7 








1.4 


00043 


7.930 


0.110 


76 ± 8 


0.21 ± 0.04 


1.2 ± 0.2 


1.7 ± 0.3 


3 


3 





3 


1.3 


00057 


7.030 


0.150 


114 ± 11 


0.21 ± 0.04 


-1.0 ± 0.2 


1.6 ± 0.3 


9 


9 








1.4 


00061 


7.680 


0.150 


93 ± 9 


0.17 ± 0.04 


-1.0 ± 0.2 


1.8 ± 0.3 


10 


10 








1.3 


00079 


7.960 


0.250 


71 + 7 


0.23 ± 0.05 


1.1 ± 0.2 


1.3 ± 0.3 


11 


11 





11 


1.3 


00103 


7.660 


0.150 


87 ± 9 


0.20 + 0.045 


-1.0 ± 0.2 


1.6 ± 0.3 


9 


9 








1.5 


00119 


8.420 


0.150 


65 ± 7 


0.18 ± 0.04 


1.2 ± 0.2 


1.9 ± 0.4 


10 


10 





10 


1.6 


00151 


9.100 


0.150 


43 ± 4 


0.22 ± 0.04 


1.1 ± 0.2 


1.4 ± 0.3 


20 


21 





21 


1.7 


00158 


9.270 


0.150 


44 + 4 


0.18 ± 0.04 


1.3 ± 0.3 


1.5 ± 0.3 


15 


15 


15 


15 


1.7 


00192 


7.130 


0.030 


99 ± 10 


0.25 ± 0.05 


1.2 ± 0.2 


1.5 ± 0.3 


8 


8 





8 


1.6 


00192 


7.130 


0.030 


99 ± 10 


0.25 ± 0.05 


1.2 ± 0.2 


1.5 ± 0.3 


6 


6 





5 


1.6 


00245 


7.820 


0.150 


80 ± 8 


0.21 ± 0.04 


1.0 ± 0.2 


1.6 ± 0.3 


7 


7 





7 


1.3 


00288 


9.840 


0.150 


33 + 3 


0.19 ± 0.04 


0.9 ± 0.2 


1.4 ± 0.3 


9 


10 


11 


11 


1.3 


00371 


8.720 


0.150 


59 ± 6 


0.17 ± 0.03 


1.1 ± 0.2 


1.7 ± 0.3 


9 


10 





10 


1.3 


00532 


5.810 


0.260 


193 ± 19 


0.22 ± 0.05 


-1.0 ± 0.2 


1.1 ± 0.2 


12 


12 








1.3 


00584 


8.710 


0.240 


54 ± 5 


0.20 ± 0.04 


1.0 ± 0.2 


3.2 ± 0.7 


5 


4 





6 


1.6 


00631 


8.700 


0.150 


52 ± 5 


0.22 ± 0.04 


1.1 ± 0.2 


1.4 ± 0.3 


21 


21 





21 


1.4 


00699 


11.720 


0.150 


13 ± 1 


0.23 ± 0.05 


1.0 ± 0.2 


2.0 ± 0.4 


4 





8 


8 


1.2 


00793 


10.260 


0.150 


29 + 3 


0.16 ± 0.03 


1.1 ± 0.2 


1.7 ± 0.3 


9 


9 


9 


9 


1.1 


00793 


10.260 


0.150 


29 ±3 


0.16 ± 0.03 


1.1 ± 0.2 


1.7 ± 0.3 


13 


13 


13 


13 


1.1 


00847 


10.290 


0.150 


29 ±3 


0.16 ± 0.03 


1.0 ± 0.2 


1.6 ±0.3 


15 


15 


15 


15 


1.4 


01036 


9.450 


0.300 


36 ±4 


0.22 ± 0.04 


1.0 ±0.2 


1.2 ±0.3 


9 


9 


9 


9 


1.3 


01036 


9.450 


0.300 


38 ±4 


0.20 ± 0.04 


1. ±0.2 


1.3 ± 0.3 


13 


13 


13 


13 


1.3 


01866 


12.400 


0.150 


8+ 1 


0.29 ± 0.06 


1.1 ± 0.2 


1.3 ± 0.3 


7 





10 


10 


1.6 


01980 


13.920 


0.150 


5 ± 1 


0.16 ±0.03 


1.4 ±0.3 


2.7 ± 0.5 


5 





7 


6 


1.7 


01980 


13.920 


0.150 


6± 1 


0.15 ± 0.03 


1.4 ± 0.3 


2.1 ± 0.4 


9 


9 


10 


10 


1.7 


11500 


18.400 


0.150 


0.8 ±0.1 


0.12 ±0.02 


1.6 ±0.3 


1.7 ±0.3 


10 


20 


20 


20 


1.4 
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Table 4. Best-fitting values of physical parameters determined for the Pallas Collisional Family asteroids 
excluding (2) Pallas with WISE observations. Negative values of tj and or R p indicate that the parameter was 
not free but fixed to the corresponding positive value. Errorbars shown are minimum estimates and correspond 
to 10% relative error for D, and 20% for r], p v and R p . 



Designation 


H 


G 


D [km] 


Pv 


n 




Rp 


Wl 


W2 


W3 


W4 


00531 


12.00 


0.15 


16 ±2 


0.10 + 0.02 


0.9 ± 


0.2 


-1.5 ±0.3 





8 


12 


12 


03579 


13.60 


0.15 


7± 1 


0.12 + 0.02 


0.9 ± 


0.2 


0.4 ±0.1 


13 


14 


13 


13 


03579 


13.60 


0.15 


7± 1 


0.13 + 0.03 


1.0 ± 


0.2 


-1.5 ± 0.3 








9 


9 


05222 


11.30 


0.15 


19 + 2 


0.15 + 0.03 


1.0 ± 


0.2 


0.4 ±0.1 


13 


13 


13 


13 


05234 


11.90 


0.15 


14 ± 1 


0.15 + 0.03 


1.0 ± 


0.2 


-1.5 ±0.3 








6 


6 


05234 


11.90 


0.15 


15 ±2 


0.13 + 0.03 


0.9 ± 


0.2 


-1.5 ±0.3 








10 


10 


05330 


12.10 


0.15 


15 ±2 


0.11+0.02 


0.9 ± 


0.2 


0.6 ±0.1 


10 


16 


16 


15 


08009 


13.70 


0.15 


6± 1 


0.14 + 0.03 


1.3 ± 


0.2 


-1.5 ±0.3 








11 


9 


11064 


12.60 


0.15 


9± 1 


0.21 ± 0.04 


1.2 ± 


0.2 


0.4 ±0.1 


16 


17 


17 


17 


12377 


12.60 


0.15 


10 ± 1 


0.15 + 0.03 


0.9 ± 


0.2 


-1.5 ±0.3 








9 


9 


14916 


13.50 


0.15 


8± 1 


0.12 + 0.02 


0.8 ± 


0.2 


-1.5 ±0.3 








8 


9 


15834 


13.20 


0.15 


9± 1 


0.11+0.02 


0.8 ± 


0.2 


-1.5 ±0.3 








9 


9 


23830 


13.50 


0.15 


10 ± 1 


0.07 ± 0.01 


1.0 ± 


0.2 


-1.5 ±0.3 





5 


12 


12 


24793 


13.80 


0.15 


7± 1 


0.11+0.02 


1.2 ± 


0.2 


-1.5 ±0.3 








12 


12 


24793 


13.80 


0.15 


7 ± 1 


0.13 + 0.03 


1.1 ± 


0.2 


-1.5 ±0.3 





10 


13 


13 


25853 


13.30 


0.15 


8± 1 


0.14 + 0.03 


1.0 ± 


0.2 


-1.5 ±0.3 








6 


6 


33166 


12.90 


0.15 


10 ± 1 


0.12 + 0.02 


1.3 ± 


0.2 


-1.5 ±0.3 








14 


14 


33750 


12.50 


0.15 


12 ± 1 


0.12 + 0.02 


1.0 ± 


0.2 


0.6 ±0.1 


14 


14 


14 


14 


36273 


12.80 


0.15 


9± 1 


0.18 + 0.04 


1.0 ± 


0.2 


-1.5 ±0.3 





16 


17 


17 


39646 


13.50 


0.15 


4.4 ± 0.4 


0.36 ± 0.07 


0.8 ± 


0.2 


-1.5 ±0.3 








10 


12 


40101 


14.10 


0.15 


6± 1 


0.10 + 0.02 


1.4 ± 


0.2 


-1.5 ±0.3 








6 


3 


44232 


13.10 


0.15 


9± 1 


0.12 + 0.02 


1.1 ± 


0.2 


-1.5 ±0.3 





14 


14 


14 


46037 


13.70 


0.15 


6± 1 


0.15 + 0.03 


1.1 ± 


0.2 


-1.5 ±0.3 








8 


8 


52229 


13.60 


0.15 


8± 1 


0.09 ± 0.02 


1.0 ± 


0.2 


-1.5 ±0.3 





11 


11 


11 


57050 


13.50 


0.15 


7± 1 


0.17 + 0.03 


0.8 ± 


0.2 


-1.5 ±0.3 








6 


4 


66714 


14.30 


0.15 


6± 1 


0.09 ± 0.02 


1.1 ± 


0.2 


-1.5 ±0.3 





7 


7 


7 


66714 


14.30 


0.15 


5.5 ± 0.5 


0.11+0.02 


1.2 ± 


0.2 


-1.5 ±0.3 








7 


7 


66803 


12.50 


0.15 


8± 1 


0.32 ± 0.06 


1.0 ± 


0.2 


-1.5 ±0.3 








15 


15 


67370 


13.70 


0.15 


6± 1 


0.15 + 0.03 


1.0 ± 


0.2 


-1.5 ±0.3 








11 


12 


67779 


12.60 


0.15 


10 ± 1 


0.16 + 0.03 


0.9 ± 


0.2 


0.4 ±0.1 


11 


13 


13 


13 


69371 


13.70 


0.15 


7± 1 


0.11 ±0.02 


1.0 ± 


0.2 


-1.5 ±0.3 








21 


20 


69931 


13.90 


0.15 


7 ± 1 


0.11 ± 0.02 


0.9 ± 


0.2 


-1.5 ± 0.3 





12 


12 


12 


82899 


13.60 


0.15 


7 ± 1 


0.13 ±0.03 


1.1 ± 


0.2 


-1.5 ±0.3 








11 


11 


87006 


13.90 


0.15 


6± 1 


0.14 ±0.03 


1.0 ± 


0.2 


-1.5 ±0.3 





7 


9 


9 


90368 


13.40 


0.15 


7± 1 


0.14 ±0.03 


1.0 ± 


0.2 


-1.5 ±0.3 








8 


8 


A0590 


14.60 


0.15 


5.1 ±0.5 


0.10 ±0.02 


1.2 ± 


0.2 


-1.5 ±0.3 








20 


18 


A 1283 


13.90 


0.15 


6± 1 


0.15 ±0.03 


0.9 ± 


0.2 


-1.5 ±0.3 








7 


7 


A3779 


14.40 


0.15 


3.9 ± 0.4 


0.21 ± 0.04 


-1.0 ± 


0.2 


-1.5 ±0.3 








8 





A9640 


13.70 


0.15 


8± 1 


0.10 ±0.02 


1.1 ± 


0.2 


0.6 ±0.1 


17 


21 


21 


21 


B3770 


14.10 


0.15 


6+1 


0.11 ±0.02 


1.0 ± 


0.2 


-1.5 ±0.3 








11 


11 
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Table 4. continued. 



Designation 


H 


G 


D [km] 


Pv 






R P 


Wl 


W2 


W3 


W4 


B8223 


14.40 


0.15 


4.7 ± 0.5 


0.14 ± 0.03 


1.4 ± 0.2 


-1 


.5 ± 0.3 








13 


12 


C3349 


14.70 


0.15 


4.1 ± 0.4 


0.14 ± 0.03 


1.2 ± 0.2 


-1 


.5 ± 0.3 





10 


14 


14 


D6038 


15.00 


0.15 


5.0 ± 0.5 


0.07 ± 0.01 


-1.0 ± 0.2 


-1 


.5 ± 0.3 








5 





D8406 


14.80 


0.15 


4.6 ± 0.5 


0.10 ± 0.02 


0.9 ± 0.2 


-1 


.5 ± 0.3 








9 


8 


E5861 


14.60 


0.15 


4.5 ± 0.5 


0.13 ± 0.03 


-1.0 ± 0.2 


-1 


.5 ± 0.3 








3 





F7914 


14.10 


0.15 


7 ± 1 


0.09 ± 0.02 


1.2 ± 0.2 


-1 


.5 ± 0.3 








7 


3 


H6413 


14.70 


0.15 


4.4 ± 0.4 


0.12 ±0.02 


1.7 ±0.2 


-1 


.5 ± 0.3 








7 


5 


18324 


14.50 


0.15 


6± 1 


0.09 ± 0.02 


1.3 ±0.2 


-1 


.5 ± 0.3 








8 


4 


K6956 


14.30 


0.15 


4.8 ± 0.5 


0.15 ±0.03 


1.3 ±0.2 


-1 


.5 ± 0.3 








4 


3 


N4076 


14.20 


0.15 


4.7 ± 0.5 


0.17 ±0.03 


-1.0 ±0.2 


-1 


.5 ± 0.3 








4 
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Appendix A: Thermal modelling of 
WISE asteroid data 



Delbo et al. 



Our aim is to model the observed asteroid flux 
as a function of a number of physical parameters 
and derive the set of parameter values that most 
closely reproduce the actually measured fluxes. 
In this work we follow the method described by 
(1201 lbl) . The set of wavelengths 



Mainzer et al 



covered by WISE (specified in Sect. [2]i allow 
us to derive up to three parameters by fitting a 
thermal model to asteroid WISE data: asteroid 
effective diameter, beaming parameter and re- 
flectance at 3.4 /mi (defined below). Within the 
wavelength range covered, the observed aster- 
oid flux consists of two components: 



77( m ) r 



(A.l) 



The thermal flux component (/ t h,,i) is the main 
contribution to W3 and W4 whereas the re- 
flected sunlight component (r Sj/ t) dominates 
in band Wl. In general, W2 will have non- 
negligi ble contributio ns from both components 
( Mainzer etallbo lib ). 

The computation of is based on the 
Ne ar Earth Asteroid Thermal Model (N EATM ; 
Harris! 1998 ; 



see 



s2 



Delbo & Harris 2002). The as 



teroid is assumed to be spherical and its sur- 
face is divided into triangular facets which con- 
tribute to the total thermal flux observed by 
WISE in accordance to the facet temperature 
(Ji), the geocentric distance (A) and the phase 
angle (a Q ). In turn, the temperature of each facet 
depends on the asteroid heliocentric distance 
(r Q ) and its orientation with respect to the di- 
rection towards the sun. It is given by 



-f-(l - A)jjjScii = TjcreTfdaj, 



(A.2) 



which results from assuming that each surface 
element da/ is in instantaneous equilibrium with 
solar radiation. S is the solar power at a dis- 
tance of 1 AU, A is the bolometric Bond albedo, 
e is the emissivity (usually taken to be 0.9; see 



2007, and references therein), cr is 
the Stefan-Boltzmann constant and /i,- = cos 
where 6>, is the angle between the normal to 
the surface element i and the direction towards 
the sun. Non-illuminated facets will be instan- 
taneously in equilibrium with the very low tem- 
peratures of the surroundings (~0K) and thus 
their contribution to ffoj is neglected in the 
NEATM. Finally, the beaming parameter (77) 
can be thought of as a normalisation or calibra- 
tion factor that accounts for the different effects 
that would change the apparent day-side tem- 
perature distribution of the asteroid compared to 
that of a per fectly smooth, non-rotating sphere 
(lHarrislll998l) . These include, for example, the 
enhanced sunward thermal emission due to sur- 
face roughness (77 < 1), or the non-negligible 
night-side emission of surfaces with high ther- 
mal inertia which, in order to conserve energy, 
causes the day-side temperature to be lower than 
that compared to the ideal case with zero ther- 
mal inertia (77 > 1). 

The asteroid thermal flux component is then 
given by 

/ th ,^ = 0^/ y (r i ), (A.3) 

i 

where f,j is the contribution from each illumi- 
nated facet of a 1-km sphere; Q = (D/lkm) 2 
scales the cross-section of the latter to the cor- 
responding value of an asteroid of diameter D. 
The color correction associated with each value 
of Ti and each WISE band is applied to the facet 
flux. By definition, it is the quotient of the in- 
band flux of the black-b ody at the given te m- 



perature to that of Vega (IWright et al 



2010). A 



color correction table was generated for all in- 
teger temperatures from 70 K up to 1000 K us- 
ing the filter profiles available from ICutri et al. 
J2OIJ). 

The reflected light component, the second 
term in the right hand side of Eq. IA.1I is calcu- 
lated as follows. First, the asteroid visible mag- 
nitude (V) that would be observed at a given ge- 
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ometry (r , A and a e ) can be est imated using We use the Lev enberg-Marquardt algorithm 



the IA U phase curve correction (IBowell et al 



19891) . along with the tabulated values of aster- 



oid absolute magnitude (H) and slope parame- 
ter (G) from the Minor Planet Center. Secondly, 
knowledge of the solar visible magnitude and 
flux at 0.55 //m (V and fy a , respectively) al- 
lows us to calculate the sunlight reflected from 
the asteroid at that particular wavelength: 

V-Vq 



r v = f Vo x 10 « 



(A.4) 



If we assume that the sun is well approximated 
by a black body emitter at the solar effective 
temperature (T = 5778 K), the estimated re- 
flected flux at any other desired wavelength (r^) 
can be computed by normalising the black body 
emission B^(Tq) to verify ry, i.e. 
Bi(T B ) 



rx = r v - 



Bv(T Q ) 

In this approximation, we can also consider 
Bir(T q ) ^ fm e 
Bv(T B ) fv 
from which we arrive at the following expres- 
sion: 



(A.6) 



(A.7) 



B, R (T Q ) 

where the subscript IR denotes 3.4 jum. Note 
that we do not color correct this component 
given the small co rrection to the flux o f a G2V 
star (see Table 1 of Wright etaklboiO ). Finally, 



in order to account for possible differences in 
the reflectivity at wavelengths longward of 0.55 
yum, a prefactor to r\ is included in the model, 
such that 



R P n. 



(A.8) 



This prefactor, R p , is by definition equivalent to 
the ratio of pm and the the visible geometric 
albedo, so we will refer to it as the "albedo ra- 
tio". The paremeter pj r is the reflectivity at 3 .4 
and 4.6 pm defined by 



Mainzer et al 



(2011b) 



To sum up, the observed model flux can then 
be written as: 



(A.9) 



dPress et al.ll 19861) in order to find the values of 
asteroid size (D = yfo., in km), beaming param- 
eter (rf) and albedo ratio (R p ) that minimise the 
X 2 of the asteroid's WISE data set, namely 



X' 



-Z 



(A. 10) 



where Fjj and crj A are the measured fluxes and 
corresponding uncertainties, /' runs over the ob- 
servation epochs and A labels the WISE bands. 
The implementation of this technique involves 
the calculation of the partial derivatives of 
with respect to the fitting parameters, which is 
straightforward in the case of Q and R p . The 
partial derivative with respect to 77 can be de- 
rived from: 



(A.5) df 



(in) 



( 8F 



(m) \ 



8Ti 



dTj 
drj 



(AM) 



The derivation of the exact analytical expresion 
for the first factor is considerably long to be in- 
cluded here, whereas the derivative of with 
respect to r\ is easily obtained from Eq. IA.2I 

Appendix B: Comparison with 



Masiero et al. 



(2011) 



Figure [2] shows that our parameter dete rmi- 
Masiero et alJ (i201ll) are 



nations and those of 



compatible in spite of the slight differences in 
the data set and the thermal modeling used in 
this work (refer to Sect. |2] and Appendix |A}, 
from which we do not expect to obtain exactly 
the same best-fit parameters for each object. In 
order to carry out a detailed c omparison be 
tween our results and those of Masiero et al. 



(201 1), we computed the mean fractional differ- 
ence (e) and corresponding standard deviations 
of D, T], pv and R p . Let e = 100(x - xm)/x, 
where x is the parameter value for a given object 
in this work and xu is t he correspoding valu e 



taken from Table 1 by 



Masiero et al. 



(2011). 



The distributions of s values are plotted in Fig. 
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Fig. B.l. Fractional difference histograms of D, rj, pv and R p . We define s = I00(x - 
wher e x is the parameter value in this work and xm is the correspoding value taken from Table 
1 by 



Masiero et al 



(1201 II) . The vertical lines mark the corresponding average values. Note: only 



parameters resulting from the same input values of H contribute to these histograms. 



IB. II These histograms only include parameter 
determinations that have the same H as input 
in order to identify possible discrepancies in re- 
sults not caused by different values of H. We 
find that our values of D and rj tend to be slightly 
greater by 1% and 3%, respectively, whereas 
our pv values are lower by 2%, though these 
deviations are small compared to the errorbars. 
On the other hand, there is a large bias towards 
lower values of R p that, while still being within 
the errorbar, must be addressed. 

Most probably, the R p discrepancy is asso- 
ciated with how the reflected flux r% is calcu- 
lated. In particular, we take the solar flux at 3.4 
pm (Jm e in Eq. IA.7l ) from the solar power spec- 
trum at zero air mass of WehrliQ, based on that 



by iNeckel & Labsl (|1984). Any differences in 
input, including solar visible magnitude, taken 
from tabulated data sources that may cause our 



r A to be systematical ly 10% greater than that of 



Masiero et al. 



(1201 II) would explain our higher 
values of R p . For instance, taking into account 
that there is only one optimum value of r s j to fit 
a given Wl data set, from Eq. IA.8l it is clear that 
a larger o will have associated a lower best-fit 
value of R p . 

The Monte Carlo estimations by the 
NEOWISE team show that the errorbars asso- 
ciated to the fitting of the data are always small 
compared to the errors inherent to the thermal 
model itself. The relative errors in diameters de- 
rived from the NEATM have been characterised 



to be ~10%-15% (Harris 



2006). From these 



http://iTedc.nrel.gov/solar/spectra/amO/wehrli 1985.new.html 



facts and the widths of the e-value distributions 
of Fig. lBTl we consider it safe to assume a min- 
imum relative error of 10% in diameter and 20% 
in beaming parameter, pv and R p . On the other 
hand, large uncertainties in the absolute magni- 
tude (sometimes as large as ~ 0.3 magnitudes) 
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will also affect the values of py, so 20% is prob- 
ably an optimistic assumption in some cases. 

Finally, we also evaluate how differences in 
the values of H result in different values of py 
and R p . We downloaded the MPC orbital ele- 
ment file as of May 2012 and compared the val- 
ues of_absolute_magnitude {Hu) to those used 
by 



Masiero et al 



(120111) . H M . About 50000 H- 
values have been updated between these two 
works, and ~38000 have been enlarged. Figure 
Q] shows a histogram of AH = Hu - Hm for the 
B-types in this work. Out of the 52 objects with 
AH + 0, as many as 43 of them have AH > 0. 
Our size determinations agree to within 10%, 
therefore larger updated values of H will result 
in lower values of geometric albedos. 



In Fig. IB. 21 we show a plot of AR p = 
R p - (R p )m versus AH for all the B-types with 
determined values of R p . The notation (R p )m 
refers to the corresp onding albedo ratios by 
(2011). There are three features 



Masiero et al. 



to note in this plot: (1) our values of R p tend 
to be ~10% systematically lower, as we already 
noted (see Fig. IB. II ); (2) most points off the 
AH = axis show a direct correlation between 
AR p and AH, as expected from the discussion 
above; (3) some points show AR p < -0.5 even 
though AH = 0. The points of feature (3) are 
ex plained by an i ncon sistency in the py values 
of 



Masiero et al 



(2011) with their correspond- 
ing values of D and H: they do not verify Eq. Q] 
and are always lower than the predicted py. 

To sum up, we have shown that if the input 
values of H are equal, our model fits are consis- 
tent within the mod el errorbars with thos e pre- 
sented in Table 1 of ] 



Masiero et al 



(120111) . The 

tendency to 10%-lower values of R p is likely 
caused by differences in solar power spectra 
data taken to estimate the reflected light compo- 
nent at NIR wavelengths (see Eq. IA.7l) . We have 
also examined how updated input values of H 
affect the best-fit parameter values and showed 
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Fig. B.2. Differences in albedo ratio determi- 
nations versus difference in absolute magnitude 
correspo nding to the B - types in this paper and 
those by 



Masiero et al 



d2011l) . Note that while 
the majority of cases verify AH = 0, our R p 
values tend to be lower (1). On the other hand, 
some of the differences in albedo ratio (2) are 
explained by the differences in the updated val- 
ues of H used in this paper. Finally, some points 
show very large differences, AR p < -0.5, in 
the albedo ratio value (3) despite the fact that 
AH = (more details in the text). 

how increasing the value of H results in greater 
values of R n and vice versa. 
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